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We formulate the Born approximation for finding resonance poles in the complex plane 
for potential scattering problems. Using the method, we study the distribution of resonance 
poles for several scattering potentials. In particular, we find for an exponential potential with 
a cutoff that the cutoff generates an infinite series of extra resonance poles below and along 
the real axis, which would not exist without the cutoff. We also find for a Gaussian potential 
that the series of resonance poles approach the imaginary axis of the complex energy plane 
from left. In other words, the real parts of the resonant eigenenergis are all negative. 



Qh| §1. Introduction 

> ■ 

^ . Resonance is a phenomenon prevalent all over physics. ^^"^^ Quantum resonance 

^ I in particular was studied intensively perhaps first in the field of nuclear physics^^"^'^^ 

, and maybe later in atomic and molecular physics. ^^^"^^^ In recent years, quantum 

resonance is attracting much attention in condensed-matter physics as well,'^^^"^''-* 
thanks to the development of nanotechnology. 

Many studies on quantum resonance focus on resonant states with long lifetimes, 
. or in other words, resonance poles close to the real axis in the complex energy plane. 

I Resonant states far from the real axis tend to be ignored on the ground that they 

are very unstable and hardly observable in experiments. 

A recent study, ^^^'^■^^ however, revealed that resonance poles far from the real 
I axis do affect resonance peaks on the real axis in the form of the Fano asymme- 

■ try. A conventional (and perhaps somewhat phenomenological) understanding of 

the Fano asymmetry is that we have an asymmetric resonance peak when a bound 
state and an energy continuum are coupled. ^^-^ In contrast, the above-mentioned 
study^^-*'^^^ microscopically expanded the transmission amplitude with respect to 
resonant states and showed that quantum interference among resonant states and 
bound states results in several types of asymmetry in resonance peaks. In particular, 
when we have two neighboring resonance poles, one close to the real axis (with a 
long lifetime and a narrow peak) and one far from the real axis (with a short lifetime 
and a broad peak), interference between these poles make the peak of the former 
resonance strongly asymmetric. This opens up the possibility of experimentally ob- 
serving resonance poles far from the real axis. 

We thereby emphasize in the present paper the importance of studying the 
distribution of the resonance poles over the entire complex plane. There have been 
quite a few studies that obtained the distribution of the resonance poles mainly for 
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exactly solvable problems of potential scattering. ^^^"^'^^ In contrast, we will here 
develop a Born approximation for finding resonance poles in the entire complex 
plane. The approximation is applicable to a wide class of solvable and unsolvable 
scattering potentials. 

Let us briefly summarize the results of the Born approximation here. We will 
show the following in the present paper: 

(i) For exponential potentials, all resonance poles are on the negative imaginary 
axis of the complex wave-number plane, or on the negative real axis of the 
second Riemann sheet of the complex energy plane. 

(ii) For exponential potentials with oscillatory modulations, the resonance poles 
which would be on the negative imaginary axis of the complex wave-number 
plane without oscillation, now split into two and shift sideways into the third 
and fourth quadrants of the complex wave-number plane. 

(iii) For exponential potentials with a discontinuous cutoff, an infinite series of extra 
resonance poles appear below and along the real axis of the complex wave- 
number plane. It is therefore dangerous to introduce a cutoff carelessly perhaps 
for numerical calculations. 

(iv) For Gaussian potentials, a series of the resonance poles approach to the line of 
argument — 7r/4 in the complex wave-number plane, or to the negative imaginary 
axis in the complex energy plane. 

We believe that the third and fourth results are of particular importance. We will 
then in the final section come up with two conjectures on the distribution of the 
resonance poles for a wider class of scattering potentials. 

For readers' benefit, we in Sec. [2] summarize the basic knowledge on the resonant 
state. We then introduce the Born approximation for resonant states in Sec. [3l Each 
of Sees. HHZl presents the afore-mentioned results (i)-(iv) in this order, respectively. 
In the final section, we summarize the paper and present two conjectures on the 
distribution of the resonance poles. We add two appendices, one to summarize the 
formulation of the Jost function in one dimension, the other to review the Jost- 
function method, a numerical method of finding resonance poles. 

§2. Definition of resonant states 

Before going into the points of the present paper, we here in the present section 
give a brief review of the definition of resonant states. Resonant states can be 
described statically or dynamically. A dynamical view of a resonant state may be 
as follows; a wave packet with a distribution of energy comes into a potential and 
gets captured in a quasi-bound state for a long time before getting scattered away 
eventually. We do not take this dynamical view in the present paper at all. 

A static view of a resonant state taken in many textbooks may be a pole of the 
S-matrix in the complex energy plane. In fact, we can equivalently define a resonant 
state as an eigenstate of the stationary Schrodinger equation 

d? 

--r^^{k; x) + V{x)(l){k; x) = k^^{k; x), (2-1) 
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Table I. Summary of definitions of resonant states, anti-resonant states, bound states, and anti- 
bound states. The eigenenergy Eras is given by the eigen- wave- number fcrcs as Eras = krea^- 



Solutions of Eq. (|2^ 


Re fcros 


Im A^rcs 


Re Ercs 


Im Eras 


Riemann sheets 


Resonant state 


positive 


negative 


any 


negative 


second 


Anti-resonant state 


negative 


negative 


any 


positive 


second 


Bound state 


zero 


positive 


negative 


zero 


first 


Anti-bound state 


zero 


negative 


negative 


zero 


second 



where we put /2m to unity. The eigenenergy E of the state is given by the eigen- 
wave- number k as E = k'^. 

The definition of a resonant state as an eigenstate is achieved under the boundary 
condition that only out-going waves exist, ^^^^'^^^ which is often called the Siegert 
condition. Suppose that we have a scattering eigenfunction in the form 



, , A(/fe)e*'=^ + 5(fc)e-*^^ for x -oo, , , 

; x) = < , , ■■. (2-2) 
C{k)e for x — > -|-oo. 



Matching conditions at the scattering potential yield in the coefficients A, B and C 
the dependence on the wave number k (and hence on the energy E = k'^). In order 
to have out-going waves only, we seek zeros of the coefficient A{k) of the in-coming 
wave, 

A{kros) = 0, (2-3) 

which is indeed equivalent to seeking poles of the S'-matrix, or the transmission 
amplitude C{k)/A{k). 

We summarize here some properties of thus-defined resonant states (Table H]) ; 
see Ref. 86) for details. Solutions of Eq. (j2-3p in the fourth quadrant of the complex 
k plane are indeed called "resonant states," whereas solutions in the third quadrant 
are called "anti-resonant states." (In some context, "anti-resonance" means a dip 
rather than a peak in the energy-dependence of the transmission probability. It is 
not related to the anti-resonant state defined here.) In a resonant state, particles 
trapped in the scattering potential leak to the infinity. An anti-resonant state is a 
time-reversal state of a resonant state. They always appear in pair. In terms of 
the energy, a resonant state is in the lower half of the second Riemann sheet of the 
complex E plane, whereas an anti-resonant state is in the upper half of the same 
sheet. Therefore, the particle number in a resonant state decays exponentially in 
time, while the particle number in an anti-resonant state grows exponentially. 

Solutions of Eq. (j2-3p in the upper half of the k plane must exist on the imaginary 
axis. They are indeed bound states. To see this, go back to Eq. (|2-2p . For a solution 
of Eq. (|2-3p , the wave function (j2-2p reduces to 

(A(fcres;x) = e^'=-I^l, (2-4) 

which further reduces to 

4>{kres = if3-,x) = e-P\^\ (2-5) 
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for a pure imaginary solution k = ifi. This is indeed a bound state if /? > 0. A 
bound state is on the negative part of the real axis of the first Riemann sheet of the 
complex E plane. 

Readers may then notice that the wave function ()2-4p diverges as |x| — > oo for 
the solutions in the lower half of the k plane, namely the resonant and anti-resonant 
states. This may appear to be unphysical, but in fact indicates that the parti- 
cles trapped in the scattering potential will go away to infinity eventually. Indeed, 
we can show^^-* that the divergence is essential for the particle-number conserva- 
tion, by noting that the divergence is balanced with the exponential decay in time, 
exp(-f| Im^JresD- 

For some potentials, a resonant state in the fourth quadrant and an anti-resonant 
state in the third quadrant collide and turn into two solutions on the negative part 
of the imaginary axis of the complex k plane. Such solutions are often called "anti- 
bound states." ^'^^ In terms of the energy, an anti-bound state is on the negative part 
of the real axis of the second Riemann sheet of the complex E plane. 



We are now in position to formulate the Born approximation for finding reso- 
nant, anti-resonant and anti-bound states. We here consider quantum-mechanical 
scattering due to symmetric potentials in the infinite one-dimensional space. All 
eigenfunctions are hence classified into even and odd functions. For spherical poten- 
tials in the three-dimensional space, we choose only odd eigenfunctions because the 
radial wave function must vanish at the origin of the semi-infinite radial coordinate. 

We start with the Schrodinger equation (|2-ip . As described above, we assume 
V{—x) = V{x). In the corresponding Lippmann-Schwinger equation 



§3. Born approximation for resonant states 





the retarded Green's function in one dimension is given by 




1 



(3-2) 



and the unperturbed wave function (j)Q{k;x) is given by 




ikx 



—ikx 



(3-3) 



where 




+ 1 



-1 



for even solutions, 
for odd solutions. 



(3-4) 



In the Born approximation, we replace the perturbed wave function <j){k; x) in the 
integrand of Eq. (j3-ip with the unperturbed wave function (j)o{k; x). In the following, 
we will be interested only in the behavior in the limit x ^ oo. We hence use only 
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the part x > x' in the retarded Green's function (l3-2p . We thus use the formula 

hm ^(k; x) = e*^^ + ae-''''' + — / e^'^^^-^'Vfx') ( e^^^^' + ae"^'^^' ) dx' 



ikx f-oo 

e^kx ^ ^^-ikx ^ ^ I v{x'){l + ae-^'''''')dx'. (3-5) 



oo 



After obtaining the asymptotic form (l3-5p for a specific potential we will 

seek resonance poles by solving Eq. (j2-3p . More specifically, since we classify eigen- 
functions into even and odd functions, we symmetrize and antisymmetrize the scat- 
tering state in the form 



_ ) '^Mky^'' + {<TB{k) + C{k)) e-^^"" for x -oo, 
- \ j^^f.^^-ikx ^ ^ aC{k)) e^^^' for x +oo, ^ ' ^ 



or after changing the normalization 

4>{k;x) -- 



ae^''"'' + a x^^^^^,, ^ e^^""' for x ^ -oo, 
+ (tC(A;) 



e*'''' + for X ^ +00. 

+crC(A;) 



(3-7) 



Therefore, the Siegert condition that we have out-going waves only gives the reso- 
nance equation in the form 

^^^■•^^^ - 0. (3-8) 



res J 

This produces the resonance poles k = /cres in the complex wave-number plane and 
hence E = Ej-^s = feres ^ in the complex energy plane. 

§4. Exponential potential 

In the present section, we discuss the resonance spectrum for a potential with 
one or multiple exponential functions. We start with one exponential, 

= -Foe-'^ol^l, (4-1) 

where Vq > and kq > 0; within the Born approximation, the spectrum for the 
case of multiple exponential functions is simply a superposition of the case of one 
exponential function. 

The single exponential (j4Tp is an exactly solvable potential. The Jost solution, 
which is defined by the behavior e*'^'^' in the limit |x| — > oo, is given by^-*'^^'*^) 

f+{k; x) = exp [i{k/Ko) In {Vq/kq^)] F {1 - 2iak/Ko) J-2ife/«„ (2y/VoP^^^/^) , 

(4-2) 

where Ju{z) is the Bessel function of the first kind. The Jost function F^{k) (which 
should not be confused with the Jost solution) is defined as Eq. ()A-9p . The resonance 
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equation (I3-8P in this case is expressed as Eq. (1A-I2p . In the present case in particular, 
the resonant states come from the fohowing part of the equation: 



0=1 -^Uik^cs/Ko (2\/^/ko) for even solutions, 
^o/'^o) for odd solutions. 



(4-3) 



These equations give solutions only on the imaginary axis, namely the bound states 
and the anti-bound states. There are infinitely many anti-bound states with 



for even solutions, 
for odd solutions. 



(4-4) 



for large lA^resl- 

We are going to show below that the Born approximation for the single exponen- 
tial function (j4-ip yields a pair of anti-bound states. Although the Born approxima- 
tion does not reproduce the fact that there are infinitely many anti-bound states, it 
does reproduce the fact that there are only anti-bound states in the lower half plane, 
no pairs of resonant and anti-resonant states. This may be remarkable considering 
that the approximation is quite naive and furthermore the imaginary axis is where 
it is at its worst; the Born approximation is considered to be accurate away from the 
imaginary axis, that is, when the particle is running fast. 

Let us now start the Born approximation. The formula (j3-5p is followed by 
straightforward calculation for the potential ()4-ip . resulting 



lim (j){k; x) 



^^-ikx ^ ^ikx 



2ik \hq kq 



a 



+ 



a 



2ik kq + 2ik 



(4-5) 



We can change the normalization of the wave function (j4-5p as in Eq. (j3-7p , but we 
can do it in a more symmetric way. Considering that singularities in the lower half 
of the complex k plane give resonant states, we split the term proportional to Vq in 
Eq. (j4-5p into two and renormalize the wave function as 



lim (f){k; x) 



a 



oc ae 



ae 



ikx 



-ikx 



2ik \kq ' kq — 2ik 
1 



1 

— + 



1 

— + 



a 



1 

— + 



a 



1 + 



2ik \kq kq — 2ik 

^0 



+ e 



ikx 



2ik \kq ' kq + 2ik 
Vb / i a 



1 

— + 



2ik \kq kq + 2ik 



1 



+ 



a 



2ik \kq kq — 2ik 



+ e 



ikx 



1 



1 



+ 



a 



2ik \kq + 



(4-6) 



We thereby have a form that is symmetric with respect to k ^ —k (except for the 
factor a). 

The resonance equation (l3-8p now reduces to 



1 + 



Vo 



1 



+ 



kq kq — 2ikj., 



0. 



(4-7) 
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It is easy to see that, within the Born approximation, the resonance equation for a 
potential with multiple exponential functions 



= - ^ y^-e-'^^^l (4-8) 



is a simple extension of Eq. ()4-7p : 



3 



1 + 2^-^ —+ _ )=0. (4-9) 



It may be noteworthy that the same equation can be produced by the method 
described in Jost's paper. ^) For the potential of the form (j4-8p . Jost expanded the 
Jost solution as 

(oo 
i n=l 

for X > and obtained recursion relations for the coefficients Cj„. It turns out that 
Cjn = OiVo^). Therefore, the Born approximation is equivalent to neglecting the 
terms n > 2. As is mentioned above, the resonant states are given by the zeros of the 
Jost function defined by Eq. ()A-9p . We can show after some algebra that Eq. ()A-12p 
(with the approximation of neglecting the terms n > 2) is indeed equivalent to 
Eq. (|¥9D . 

Now we solve the resonance equation (l4-7p . We could solve it separately for 
a = +1, when the equation is quadratic, and for a = —1, when the equation is linear. 
We here solve it more systematically. We note that the dominant contribution of 
Eq. (|4-7p comes from the term with the denominator kq — 2ikres- We therefore have 

2iA;res(2iA;res - Ko) = aVo, (4-11) 

whose zeroth-order solution is 2ikj-cs — kq. After putting 2ikres = 1^0 + ^, we then 
have 

{KQ + 5)5 = aVo. (4-12) 

Since 5 = 0(Vo), Eq. (j4-12p actually reduces to hq5 = crVo in the first order of Vq. 
We thereby arrive at a pair of solutions on the negative imaginary axis: 

feres = ( + ) • (4T3) 
2 V Ko/ 



For the potential with multiple exponential functions, Eq. (j4-8p . there is a pair of 
solutions for each exponential: 

feres = —77 \ f^j + '^~] foi' each j. (4-14) 
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-0)2/2 -coi/2 



1/2 0)2/2 




Fig. 1. A schematic illustration of the resonance spectrum within the Born approximation: (a) for 
a potential with multiple exponentials, we have a pair of even and odd solutions corresponding 
to each exponential; (b) for a potential with multiple oscillatory exponentials, each solution 
splits into a resonant state and an anti-resonant state. 



This is illustrated in Fig. [I||a); we have a pair of even and odd solutions on the 
imaginary axis for each exponential in the potential (|4-8p . Note that there are no 
exact solutions for potentials with multiple exponential functions. 

We have thus showed within the Born approximation that the resonant solutions 
are only on the negative imaginary axis for potentials with one or more exponentials. 
At the end of the present section, let us address the reason why we can reproduce 
only a pair of solutions by the Born approximation. We show that the pair of 
solutions (|4-14p according to the Born approximation actually corresponds to the 
pair of the case n = 1 of the exact solutions (l4-4p . The exact solutions with n > 2 
emerge because of higher-order effects. 

For Vo ^ kq'^, the Bessel function in Eq. (j4-3p for odd solutions can be expanded 
in the form 

OC> / -1 N 77, Xl 

= U2a) = a^l' Yl (4-15) 
^-^ nil 1/ + n + 1 

n=0 

where 1^ = — 2i/cres/'*o and = Vq/kq'^. Suppose that we terminate the summation 
in Eq. (j4T5p at n = 1, that is, at the order 0(Vo), and assume the deviation of 
from the value —1 as ly = —1 + 61. The resonance equation (j4-15p then reduces to 

This is followed hy 61 = v + 1 = a, which is indeed the solution (|4T2p by the Born 
approximation. Next we terminate the summation in Eq. (j4T5p at n = 2, or at the 
order O(Vo^), and assume v = —2 + 62- The resonance equation (|4-15p then reduces 
to 



1 a 1 



a2 



(-l + ,52)-Q + 



(4-17) 



r(-i + 52) r{62) 2r{i + 62) r{62) 

We numerically confirmed that the solution of Eq. (j4-17p well reproduces the exact 
solution (l4-4p with n = 2 for a <C 1. We went further and confirmed this correspon- 
dence up to n = 5. This shows that the exact solutions (j4-4p with n > 2 are indeed 
higher-order effects. 
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Summarizing the above discussion, we stress the following point. It is not a 
coincidence that the Born approximation reproduced the resonant solutions on the 
negative imaginary axis. On the contrary, the Born approximation successfully re- 
produced the first pair of the anti-bound states. 



§5. An exponential potential with an oscillation 



In the present section, we discuss the resonance spectrum for an exponential 
potential with an oscillation; 



V{x) 



-Vo cos{uJox)e 



(5-1) 



where Vq > and kq > 0. There is no exact solution for this potential. We discuss 
this case as a primer to the case in the next section, where we discuss an exponential 
potential with a cutoff. We will argue in the next section that a cutoff plays a role 
of superposition of various oscillations. In the present section, we are going to show 
that an oscillatory factor in Eq. (j5-ip makes each of the anti-bound states (j4-13p split 
into a pair of resonant and anti-resonant states. Then in the next section, we will 
show that a cutoff, because of its multiple oscillatory Fourier components, yields a 
sequence of pairs of resonant and anti-resonant states. 

We can rewrite the potential (j5-l|) as a sum of two exponential functions;^^-* 



V{x) 



Vo 



By defining 



Kl = K0+ ilVo, K2 = HQ - ^^Q, Vl = V2 



(5-2) 



(5-3) 



we can reduce the potential (j5-2p into the form 



V{x) = -^F^-e-'^^l^l. 



Hence the resonance equation is given by 

2 



1 1 

+ cr- 



0. 



(5-4) 



(5-5) 



The solutions within the Born approximation for multiple exponentials are given 
by Eq. (j4-14p . In the present case of Eq. (|5-3p , we readily obtain the solutions as 



kres = -7; [ f^0± i^^O + 

2 \ 2 Ko± lU)o 



1 ( Vo/2 

2 "V Ko2+L^oV 2 



(5-6) 
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V{x) 




Fig. 2. An exponential potential with a cutofF. 



The answer converges to (|4-13p as loq ^ 0; note that we should also replace Vq/2 
with Vo to make the exact correspondence. The largest difference between Eq. (j4-13p 
and Eq. (I5-6P is the existence of the real part. Each of the anti-bound states given by 
Eq. (I4-13P now splits into a pair of a resonant state and an anti-resonant state and 
shifts sideways. If we have more than one exponentials with different oscillations, 
each anti-bound state splits into a pair as illustrated in Fig. mb). 



§6. An exponential potential with a cutoff 



In the present section, we discuss the resonance spectrum for an exponential 
potential with a cutoff (Fig. [2]); 

{-"Kie"""!''! for Ixl < b 
° - ' (6-1) 

for |x| > b, 

where Vq > 0, > and 6 > 0. This potential is exactly solvable; we nonetheless 
treat the problem in the present framework of the Born approximation. The purpose 
is to see explicitly the origin of a resonance spectrum surprisingly different from the 
case of a simple exponential potential. 

As we described in section HI an exponential potential produces only anti-bound 
states. In striking contrast, an exponential potential with a cutoff produces a series 
of extra resonant states below and along the real axis of the complex wave-number 
plane. Let us start with somewhat hand-waving argument, where we attribute this 
structure to an exponential with various oscillations. The discontinuity at the cutoffs 
X = ±5 may be described as superposition of various Fourier components. In other 
words, the cutoff potential (j6-ip may be written in the form 



V{x) ~ -yoe"''°l'^l ^CnCOSWnX. (6-2) 



(6-3) 



If we accept this, Eq. (j5-6p dictates that we have a series of resonant poles at 

±UJn - iHO 
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These poles are aligned below and along the real k axis. 

For a little bit more detailed analysis, let us express the potential ()6Tp as 



V{x) 



-Voe 



-K0\X\ 



X {0{x + b) - 0{x - b)) , 



(6-4) 



where 0{x) is the Heaviside step function. The Fourier expansion of the term with 
the step functions is given by 



0{x + 6) - 0{x - b) 



sin kb 



(6-5) 



The Fourier coefficient sinkb/k has peaks and dips at A; = (2n — l)7r/(26) with an 
integer n except for k = ib7r/(26). Thus we could argue that the approximation (j6-2p 
is legitimate with 

2n - 1 , , 

(6-6) 



iOr, 



2b 



-vr 



for large |n|. Indeed, we will see in Eq. (j6-20p below that the location of the poles 
can be approximated as 

(2n - l)7r 



46 



(6-7) 



for large |n|. This is consistent with Eq. ()6-3p accompanied by Eq. (j6-6p . 

Let us now develop the Born approximation for the exponential potential with 
a cutoff, Eq. (j6-ip . The formula (j3-5p gives 



lim 

X—I-+00 



:x 



1 



2ik 

,-Kob 



kq kq — 2ik kq + 2ik 



2ikb 



+ 



ere 



-2ikb 



kq kq — 2ik KQ + 2ik 
After the procedure given in Eq. (l4-6p . we have the resonance equation as 



0=1 + 



1 + 



1 

— + 



ae 



1 



-Kob 



+ a- 



1 



Kq Hq 2ikj-i,g 

g-(K0-2ifcrcs)fc 



Kq — 2ikj-. 



(6-9) 



In order to solve this resonance equation, we again note that the dominant 
contribution comes from the term with the denominator kq — 2ik-res- We therefore 
have 



(JVq 



1 



j2ikias-Ko)b 



(6-10) 



For the equation to converge in the limit 6 — > oo, the real part in the exponent must 
not depend on b. We therefore assume the form 

a + i(3 



(6-11) 
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where a and (3 are real. Using the form (j6-ll|) in Eq. (j6-10p for large b, we have 



a + i 



or 



a — a- 



-a e cos p, 

Ko 



Vob a ■ o 

-a e smp. 

Ko 



Taking the ratio of the two, we obtain 

tan (3 



P 



a — aVob/Kf) 
whose solution is approximated for large by 

2n- 1 



-vr + 5n, 



(6-12) 



(6-13) 



(6-14) 



(6-15) 



where n is an integer and 6^ denotes the correction. Equation (|6-13p is then approx- 
imated by 



, Vob 
±a e dn, 



Vob 
a — a 

Ko Ko 

2n - 1 Vob 

TT ~ ±0" e 

2 Ko 



(6-16) 
(6-17) 



where, on the right-hand sides, we choose the positive signs for odd n and the negative 
signs for even n. (Since the left-hand side of Eq. (j6-17p is negative for /3 > 0, we 
realize that in the region ReA;res > 0, the symmetric solutions a = 1 appear for even 
n and the antisymmetric solutions a = —1 appear for odd n. The situation is the 
opposite for Re/cres < 0.) We then have 



a = In 



Ko |2n — 1| 
2 



-vr 



(2n - l)7r 



a- 



.Vob 

Ko 



In 



Ko |2n — 1| 

Vob : 



-TT 



(6-18) 
(6-19) 



To summarize, we arrive at the solutions 



2n- 1 
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i 
2 



-vr + 



1 



(2n- l)7r 



' Vo 1, 
a in 

Ko b 



Ko |2n — 1| 

Vob ; 



-vr 



Ko + 7 In 
b 



Ko |2n — 1| 

Vob : 



-vr 



(6-20) 



for large b and n, or, in the lowest order, Eq. (j6-7|) . As we emphasized above, this is 
consistent with the solutions for the exponential potential with oscillations generated 
by the two step functions (j6-4p . In short, the introduction of a cutoff possibly for 
numerical calculations can be disastrous, generating infinitely many extra poles that 
would not exist without the cutoff. 
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Fig. 3. Exact resonant solutions of an exponential potential with a cutoff with Vb — 10 and b = 
1,2,5, 10[/«:o]. 



At the end of the present section, we show in Fig. [3l the results of the exact 
solution for the exponential potential with a cutoff, Eq. (16-ip . and compare them 
with the result of the exact solution for the exponential potential without a cutoff, 
Eq. ()4-ip . The results are consistent with the Born approximation (j6-20p particularly 
for large b. The series of poles below and along the real k axis become denser 
and denser as we increase the cutoff range b. In fact, the residues of the poles 
decrease exponentially as the discontinuity at the cutoff exp(— kq^) decreases. Since 
the density of the poles increases only proportionally to b, the whole contributions 
of the resonant and anti-resonant poles vanish in the limit b ^ oo and hence only 
the anti-bound states remain for the exponential potential without a cutoff. This is 
also observed in the transmission coefficient shown in Fig. HI The resonance peaks 
(or dips) become less prominent as we increase b and eventually indiscernible from 
the case without a cutoff. 



§7. Gaussian potential 



In the present section, we discuss the resonance spectrum for a potential with 
one or multiple Gaussian functions. We start with a potential with one Gaussian 
function 

y(x) = -yoe-""'"', (7-1) 

where Vq > and oq > 0. We will show that the series of resonant eigen-wave- 
numbers asymptotically approaches the line Im k = — Re k from below, or the line of 
argument — 7r/4 in the complex k plane from below, and the imaginary axis Re-E = 
from left. 





Fig. 4. The transmission coefficient of the exponential potential with a cutoff with Vo = 10 and 
b — 1, 2, 5, 10[/ko] (from the right curve to the left curve in this order). The panel (b) is a 
blowup of the panel (a). The case buo = 10 is indiscernible from the case without a cutoff. 



The formula (|3-5|) for the Born approximation in this case yields 



2iK ao 



(7-2) 



Distribution of Resonant Eigenvalues of Quantum Potential Scattering 15 



By changing the normahzation, we have 



Hm (/)(/c; x) oc e*'^^ + ae 



-ikx 



a;— >+oo 



^0 VTT ^ , _„-fc2/„g2 



which is followed by the resonance equation 



ao 



1 + c7e-'=^<=='/"o' 



0. 



(7-3) 



(7-4) 



Within the Born approximation, the resonance equation for the potential with mul- 
tiple Gaussians 



Vix) 



Y^V.e-^^ 



(7-5) 



is given by 



2ik,es + ( 1 + ere" 



'^rcs / 



2 /„ .2 
3 



(7-6) 



Let us now solve the resonance equation ()7-4p for large | kj-^s \ . For the resonance 
equation (j7-4p to be stable, the real part of the exponent k^es^ /qq^ should stay small 



Re 



k 2 



0(1) 



for large | A:i.es | ■ This means that if we assume the form 

^res 

ao 



a — if3 



(7-7) 



(7-8) 



with positive a and /3, we should have a — [5 = 0(1). This is then followed by 
2ikj-(.s oc 1 + i. Equation (j7-4p then gives 



exp 



k 2 
"res 



OC — (t(1 + i) 



with a positive proportional coefficient, or 



Im 



kr, 



2a(5 ~ ( n 



(7-9) 



(7-10) 



for a large integer n. We choose an even n for a symmetric solution a = +1 and an 
odd n for an antisymmetric solution a = —1. Let us then denote the solutions of a 
and P as q„ and Pn, respectively, for the integer n. We can assume the following: 



Or, 



an ^ Pn 



-In, 



(7-11) 
(7-12) 



16 



N. Hatano 



where 7„ is a real number of order nP. We will then ask for the behavior of 7^. 
Equations ([7TT|) and ([712]) give 





/ ^res \ 






exp 




= exp 


In 











7n + n 



a 

7i 



(l + i)e^", 



where on the right-hand side we use the fact that a = +1 for even n and a 
odd n. The resonance equation (j7-4p now reduces to 



2^ - -^e^" 



0. 



(7-13) 
— 1 for 

(7-14) 



We thus have the asymptotics 



7n 



■ Inn. 



(7-15) 



Note that the asymptotics ()7-15p is compatible with the asymptotics ()7-12p . because 
In n is of lower order of n than y/n. 

We compare in Fig. [5] the results of the Born approximation with the results 
of the Jost-function method, which is described in Appendix B. The Jost-function 
method is supposedly producing numerically exact solutions here. We can see that 
the results of the Born approximation are almost equal to the numerically exact 
solutions for Vq = 0.1 and fairly consistent with them in the region Re/c/ao ^ 3 for 
large Vq = 10. 

The resonant poles approach the line of argument — 7r/4 from below in the com- 
plex k plane, which is also consistent with the result of the Born approximation, 
7„ > 0. In the case of the strong potential, Vq = 10, the resonant poles are much 
closer to the — 7r/4 line, but they still approach the line from below. This means that 
all the resonant eigenenergies have negative real parts: 

Re E^es = Re /c^es^ < 0. (7-16) 
Finally, for a positive potential (Vq < in the present notation), Eq. ()7T3p reads 



exp 



Tes 

ao 



exp 



7n + ?^ + 



1 



a 

71 



(1 + i)e'^" 



(7-17) 



with a = +1 for even n and a = —1 for odd n. In other words, the even solutions 
and odd solutions are interchanged with each other. Equation (j7T4p reads 



2^ - 11:^6^" ~ 



(7-18) 



with the same Eq. (j7-15p . The fact that the resonant poles approach the line of 
arugmet — 7r/4 from below remains the same. 
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Fig. 5. The distribution of resonance poles in the complex wave-number plane obtained by the 
Born approximation and the Jost-function method, (a) ao = 1 and Vo = 0.1 and (b) ao = 1 and 
Vo = 10. In the latter case, we have three bound states. 



§8. Summary and discussions 

In the present paper, we formulated the Born approximation for finding reso- 
nance poles in the complex plane. We found it working well for scattering potentials 
of exponential functions as well as Gaussians. 

The greatest merit of using the Born approximation is that we can see the origin 
of the structure of the resonance spectra. In the case of the exponential potential 
with the cutoff in Sec. [U we were able to see clearly that the reason why the extra 
resonance poles are aligned below and along the real k axis is because the cutoff 
generates oscillatory potentials. 

The hand- waiving argument around and below Eq. (|6-2p . which was confirmed 
by the Born approximation, should be applicable to any scattering potentials with 
discontinuities. We thereby propose the following conjecture: 

Conjecture 1: Any discontinuities in the scattering potential cause a 
series of extra resonance poles below and along the real k axis as well as 
along the real E axis. 
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We were also able to explain clearly the reason why the resonance poles approach 
the — 7r/4 line in the complex k plane from below for the Gaussian potential. The 
reason of the convergence to the — 7r/4 line in the complex k plane is the fact a — P = 
0(1), which follows from Eq. (j7-7p . We note that the square of kres in Eq. (j7-7p . or 
in the exponent of the resonance equation (j7-4p . is originated from the square of x 
in the exponent of the Gaussian potential (|7-ip . It is then reasonable to conjecture 
the following; 

Conjecture 2: For a scattering potential of the form 

V{x) = -yoe""°'l^'l\ (8-1) 

the resonance poles should approach the line of the argument — 7r/(2A) in 
the complex k plane and the line of the argument — vr/A in the complex 
E plane. 

As we saw in SecSl we indeed have poles only on the imaginary k axis for A = 1, 
namely for the exponential potential. We could also regard that the limit A — > oo 
corresponds to a potential with a very rapid decrease, or a cutoff, in which case the 
poles are indeed aligned along the real k axis. 
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Appendix A 

Jost solutions and Jost functions in one dimension 

In this section, we derive the Jost solutions and the Jost functions in one di- 
mension. In the usual context, they are derived in the radial coordinate of the 
three-dimensional space, which is a semi-infinite one-dimensional space. We here 
consider the one-dimensional space which is infinite in both positive and negative 
directions. We hence suppose that reviewing the derivation for the one-dimensional 
space may be instructive. 

When the potential is symmetric with respect to the origin, any solutions can be 
classified into even or odd solutions. In fact, the solutions in the radial coordinate 
are equivalent to odd solutions in the infinite one-dimensional space, because they 
must vanish at the origin. 

We define the Jost solutions f±{k;x) by their behavior in the far right: 

lim e^^'=^/±(fc;j;) = 1. (A-1) 

X— >oo 

For a symmetric potential, the Jost functions are also classified to even and odd 
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solutions. Therefore, the Jost solutions behave in the far left as 

hm e±^^^ f^ik; x) = a ^ {^^ (A.2) 
x^-oo 1—1 for even solutions. 

The Jost solutions are generally discontinuous or undifferentiable at the origin. 
On the other hand, a physical solution (f){k; x) must be continuous and differentiable 
at the origin. The boundary conditions for the physical solution at the origin can be 
given in the forms 

0(fe;O) = l and 0'(A;;O) = O for even solutions, (A-3) 
(^(fe;0) = and (/>'(A;;0) = 1 for odd solutions, (A-4) 

where 4>'{k; x) = (90(/c; x)/dx. 

Since the Jost solutions f±{k;x) are mutually independent solutions of a second 
differential equation, we can express the physical solution as a linear combination of 
the Jost solutions: 

4^{k;x)=C+{k)f+{k-x) + C.{k)f.{k;x). (A-5) 

In order to proceed, we use the Wronskian 

W{f,g) ^ f{x)g'{x) - g{x)f'{x), (A-6) 

which must be independent of x for linearly independent solutions / and g of the 
Schrodinger equation (12Tp . Since we can obtain W{f+, f-) = —2ik in the limit 
X — > oo, we have 

Wif+,<l>) = C-{k)W{f+J.) = -2ikC^{k), (A.7) 
= C+{k)Wif-,f+) = 2ikC+{k) (A-8) 

in the limit x — > oo. In the limit x — > 0, on the other hand, we have 

F±{k) = W{f±, cf)) = f±{k- O)0'(A;; 0) - ^ Q)f'^{k- 0) 
I —f'j^{k; 0) for even solutions, 



f±{k; 0) for odd solutions, 



(A-9) 



where we used the boundary conditions ()A-3P and ()A-4p in the second equation. The 
functions F±[k) defined above are called the Jost functions, which should not be con- 
fused with the Jost solutions f±{k; x). (The Jost function for even solutions may not 
be common because most papers discuss only radial functions in three dimensions, 
in which case the solutions must vanish at the origin. ) Equating Eqs. (TATjl - dA^ . 
we obtain C±{k) = zizF^{k)/{2ik). The physical solution (IA-5P is now given in the 
form 

= ^{F-{k)f+ik;x) - F+{k)f-{k;x)) . (A-IO) 
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The Siegert condition for a resonance pole is 

lim (/)(fcres;a:;)e"*''"''="'^' = constant. (A-11) 

a;|^oo 

In other words, we should have only the first term on the left-hand side of Eq. ()ATOp . 
Therefore, the resonance equation is given by the zeros of the Jost function, 

F+(/Crcs)=0. (A- 12) 

Appendix B 

Jost- function method of finding resonant states 

In the present section, we describe some details of the Jost-function method, 
which we used to obtain the results in Fig. [5l The method was proposed in Refs. 89)- 
91) and was extended in Ref. 23). They all obtained the solutions for Coulombic 
scattering potentials with an angular momentum. Since we do not include a Coulomb 
potential nor an angular momentum in the present paper, the formulation becomes 
simpler than the one described in the references. We consider that it may be worth 
reviewing here the simpler version. 

Let us find the solutions of the Schrodinger equation (|2-ip of a form other than 
Eq. ()A-10p . If we did not have the potential V{x), the solution would be given by a 
linear combination of the simple plane waves e^*'^^, 

^(k;x) = ^ (F^{k)e''''' - F+{k)e-'''''A fovV{x) = (B-1) 

with the Jost functions F±{]i) for V{x) = being the coefficients. When we turn 
on the scattering potential V{x\ we can write down a solution by giving the two 
coefficients the x-dependence: 

4>(k;x) = \(T-{fz;xy^'' -T^{k;x)e-'^''\ forF(x)/0. (B-2) 

That is, instead of modifying the plane waves e^*'^^ to the Jost solutions f±{k;x) 
as in Eq. (jATOp . we modify the coefficients from F±(k) to J-±{k]x). By comparing 
Eqs. (IXTO]) and (|B^ . we notice 

Mui T±{k-x) = F±{k). (B-3) 

x~*oo 

Therefore, the resonance equation ()AT2p is now converted into 

lim ^+(A;,es;x) = 0. (B-4) 

The remaining task for us is to derive differential equations for the functions 
T±{k\x). There is, in fact, redundancy in the functions ^±(/c;x), something like a 
gauge degree of freedom. In order to fix the gauge, we assume a condition 

T'„{k-x)e'^'' - T'{k-x)e-'^'' = Q, (B-5) 
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where J^^{k;x) = d!F±{k;x) / dx. 

By differentiating Eq. (|B-2p . we have 

<t^'{k-x) = ^ {j^'_{k-.x)e'^^ - :F'+{k-x)e-'^^) 

+ ^ (^-(A;; x)e'^- + T+{k- x)e-'^^) (B-6) 

= ^ {j'-{k;x)e'^- + T+{k-x)e~^^-) , (B-7) 

where we used the gauge- fixing condition ()B-5P in moving from Eq. ()B-6P to Eq. ()B-7p . 
By differentiating Eq. ()B-7p . we have 

ct>"{k;x) = ^ [T'_{k;x)e'''-+T'_,{k;x)e-'''^) 

+ ^ (.F„(A;;x)e^'=^ -.F+(fc;x)e-*'=^) (B-8) 

= ^ (.Fi(fc;2;)e^'=^+.F;(A:;x)e-^'=^) - k^ (j){k; x) . (B-9) 
Comparing this with the Schrodinger equation (j2-ip . we have 

^ (.F:(A:;x)e'^" +.F;(fc;x)e-'*=") = ^ (.F_(A:; x)e'^" - .F+(A;; x)e-'*^") . (B-10) 

Equations ()B-5P and ()B-10p constitute a set of first-order differential equations 
for the functions J^±{k;x). By solving them with respect to T^{k;x), we arrive at 

^4(A:;x) =6^*^^^^ f^_(A:;x)e*'=^-.F+(A;;x)e-*^^) . (B-11) 

The boundary conditions ()A-3p and ()A-4p are now translated into the conditions 

—J--^{k;0) = J^^(k;0) = ik for even solutions, (B-12) 
J'+ik; 0) = J^-{k; 0) = 1 for odd solutions. (B-13) 

We first solve the set of differential equations ()B-lip from x = with the bound- 
ary conditions (|B-12p and (|B-13p up until a large enough x = Xmax, scanning the 
complex k plane. We then make a density plot of |.7-+(A;; Xmax)! in the complex k 
plane. The density plot shows dips, which give rough estimates of the zeros of the 
function J^^(k; x^g^x)- Then we can use the standard Newton method to find each 
zero more precisely. In other words, we numerically solve the equation 

.F+(A;res;Xmax) =0 (B-14) 

instead of Eq. ()B-4p . This is how we obtained the results in Fig. [SI where we used 
Xres = 15/ao after confirming that the function T^(kres] Xraax) converges to a constant 
for X > 15/ao. 
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